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ABSTRACT 

We describe Spitzer IRS spectroscopic observations of the ^ 10 Myr-old star, EF Cha. Composi- 
tional modeling of the spectra from 5 /xm to 35 confirms that it is surrounded by a luminous debris 
disk with L^/L^ ^ 10~^, containing dust with temperatures between 225 K and 430 K characteristic 
of the terrestrial zone. The EF Cha spectrum shows evidence for many solid-state features, unlike 
most cold, low-luminosity debris disks but like some other 10-20 Myr-old luminous, warm debris disks 
(e.g. HD 113766A). The EF Cha debris disk is unusually rich in a species or combination of species 
whose emissivities resemble that of finely-powdered, laboratory-measured phyllosilicate species (talc, 
saponite, and smectite), which are likely produced by aqueous alteration of primordial anhydrous 
rocky materials. The dust and, by inference, the parent bodies of the debris also contain abundant 
amorphous silicates and metal sulfides, and possibly water ice. The dust's total olivine to pyroxene 
ratio of ^ 2 also provides evidence of aqueous alteration. The large mass volume of grains with sizes 
comparable to or below the radiation blow-out limit implies that planetesimals may be colliding at 
a rate high enough to yield the emitting dust but not so high as to devolatize the planetesimals via 
impact processing. Because phyllosilicates are produced by the interactions between anhydrous rock 
and warm, reactive water, EF Cha's disk is a likely signpost for water delivery to the terrestrial zone 
of a young planetary system. 

Subject headings: astrochemistry,infrared:stars, planetary systems:formation, planetary systems: pro- 
toplanetary disks, radiation mechanisms:thermal, techniques: spectroscopic, stars: 
individual (EF Cha) 



1. INTRODUCTION 

The 10-20 Myr age range is an important time in 
the formation of icy planetesimals and planets around 
^ 2 Mq stars. While debris emission decays with 
tJ" , time for t > 20 Myr as expected for simpl e collisional 
grinding of a fixed-mass planetesimal belt (|Rieke et al.l 
I 120051 : IWvatt et ani2007f ). debris emission for A/F stars 
, is observed to ris e from 5 My r to 1 Myr and peak 
from 10-20 Mvr (ICurri e et al. ' 2008al: iHernandez et al.l 
120091 : ICurrie et all l2008b. .200C The rise in debris 
emission has been attributed to the agglomerate ac- 
cretion of 1000 km-sized icy bodies in outer disk re- 
[ gions as planet formation switches from a fast runaway 
■ growth stage to a slower, debris-produ cing oligarchic 
C^ j stage (|Kenvon and BromlevI l2008l [2010D . Dynamical 
perturbations by massive planets may also explain the lu- 
minous debris emission at these a ges (jMustill and Wvat3 
[20091 : IKennedv and WvattJ [20ll). 

Spitzer has observed a number of luminous debris disks 
in this critical stage. A wide range of processes appear to 
be responsible for their extreme infrared emission, poten- 
tially revealing a variety of situations prevailing in their 
terrestrial planet zones. For example, HD 113766A, an 
F3/5 star in 10-16 Myr-old Lower Centaurus Crux, shows 
evidence for multiple debris belts, including a warm (T 
^ 440 K) silicate dominated debris disk with a compo- 
sition consistent with S-type asteroids and whose flux 
densities require an enormous reservoir (M > MMars) of 
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colliding inaterial ()Chen et al.l 120061 : iLisse et al.l 12008') . 
iLisse et all (|2008[ ) interpreted these spectral features as 
evidence of terrestrial planet-forming events. HD 172555, 
a ^12 Myr-old A5 star, also exhibits luminous warm de- 
bris disk emission at L^j/Lh. ~ 10~^, but contains spec- 
tral features consistent with silica and SiO gas produced 
by a giant hypervelocity collision, similar to that likely 
respo nsible for forming the present-day Earth-Moon sys- 
tem (jLisse et al.ll2009t i. Eta Corvi, an 1 Gyr old F2 
star, shows an extended cold Kuiper Belt disk due to 
collisions driven by dynamical stirring, and a separate 
belt of warm inner system dust due to the collision of 
a KBO with a planet in the terrestr ial habitable zone 
(|Wvatt et all 120051: [Lisse et al.ll2011al) . 

Spatially-resolved images of other luminous debris 
disks at these ages - such as HR 4796A (AOV) and 
(3 Pic (A5V) - find evidence for debris disk dynami- 
cal sculpting, such as asymmetries or war ping, plausibly 
due to the presence of massive planets (Mouillet et alj 
19971: iGoli mowski et all [20061: Schneider et al. 200^^ 
Lagrange et al.ll2010l ). However, unlike HD 113766 A and 
HD 172555, they have rather featurele ss IR spectr a, lack- 
ing evidence for massive collisions (jChen et al.l [2006) . 
The disk around 49 Ceti (AOV) extends at high surface 
brightness to a radius of 900 AU but appears to have been 
cleared within 20 AU of the star , again possibly due to 
an unseen planet (IWahhai et al.l i2007'). Modeling of the 
images suggests that the mid-infrared emission from the 
49 Ceti di sk is dominated by s mall, short-lived grains 
0.01 Atm) (jWahhai et al.l[200l . 

EF Chameleonis (hereafter EF Cha) provides another 
interesting example of the variety of disk properties at 
this stage. It is an 10 Myr-old A9IV/V star 95 pc 
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from the Sun identified from the IRAS sateUite as hav- 
ing infrared excess emission consistent with a luminous 
debris disk of Ld/U ~ IQ-^ (for L*/ Lq ~ 9.7). A 
combination of MSX, ground based, and Spitzer data 
show that its disk contains warm dust at temperatures 
~ 240 K (jRhee et al . 2007). It is among the relatively 
rare disks whose emission is dominated by warm dust 
and, therefore, by material in the terre stri al-planet zone 
(iRhee et all [20071: iCurrie et all l2007allbl : [Morales et all 
[2009^! 

In this paper, we present and analyze high signal-to- 
noise 5-35 /um spectra of the EF Cha debris disk ob- 
tained with Spitzer. Section 2 describes our observations 
and data reduction. In Section 3, we compare EF Cha 
to a simple model for blackbody disk emission and con- 
firm that much of its IR excess is most plausibly due 
to debris emission originating from EF Cha's terrestrial 
planet-forming region. We then compare the EF Cha 
spectrum to spectra from other debris disks, protoplane- 
tary disks, comets, and sophisticated mineralogical mod- 
els in Section 4 to determine the chemical composition 
of the EF Cha debris disk dust and constrain the likely 
origin of its material. We conclude by placing the spec- 
tral analysis of the EF Cha debris disk in the context of 
mineralogical studies of the early solar system and dis- 
cuss processes likely responsible for explaining the EF 
Cha disk's spectrum. 

2. OBSERVATIONS AND DATA REDUCTION 

EF Cha was observed by the Spitzer Infrared Spectro- 
graph (IRS) on September 14 and October 12, 2008 in 
staring mode (PID 50150; AORs r25679360, r25678080, 
and r25678336). The low- resolution spectra covered 5- 
38 ^m, using all four low-resolution orders (Short-Low 

1 & 2, Long-Low 1 & 2), with exposures consisting of 

2 cycles with a ramp duration of 6 seconds each. The 
high-resolution spectra covered 9.9-19.6 iim in two cy- 
cles with a ramp duration of 30 seconds each. For the 
high-resolution spectra only, we obtained separate back- 
ground observations. 

We processed all spectra using the Spectroscopic 
ModeUng An alysis and Reduction Tool (SMART, 
iHigdon et al.|[2004l) . For the low- resolution spectra, we 
first median combined all frames for a specific setting 
(e.g. Short-Low 1) at a given nod position. Then, we 
subtracted the median combined frame for one nod posi- 
tion from each frame corresponding to the other nod po- 
sition, producing skysubtracted and bad pixel mitigated 
focal plane images. We then grouped the frames together 
by calibration file and extracted spectra from them simul- 
taneously, using the automatic optimal extraction algo- 
rithm provided by SMART. Finally, we grouped the ex- 
tracted spectra for each setting together, masking the 
bonus order, clipping obvious outlying pixels and others 
using the SMART "special" clipping algorithm set at a 
2.5 (7 threshold, masking pixels lying outside the wave- 
length range where IRS calibration is reliable (e.g. A > 
38 fJ-m), and median combined all spectra. 

For the high-resolution spectra, we first median- 
combined the science frames and the sky background 
frames at a given nod position. The appropriate median- 
combined sky frame was then subtracted from the sci- 
ence frames at each nod position. Using the full aper- 
ture algorithm, we extracted each sky-subtracted frame. 



The spectra were then grouped together; obvious outly- 
ing pixels are identified by comparing pixel values from 
the four separate exposures, masked, and removed. Each 
spectrum was defringed with the algorithm irsfringe 
ported into SMART and trimmed of pixels with unreli- 
able flux calibration as determined from Table 5.1 in the 
IRS data handbook. As with the low-resolution spectra, 
we clipped the processed spectra at 2.5 a and median 
combined them to obtain the final high-res spectrum. 

Figure [T] (top panel) shows the extracted spectra as 
a function of wavelength. The low-resolution spectrum 
(black line) appears to be smooth and high quality from 
5 through ^ 30-35 fim. As expected, the high-resolution 
spectrum (grey line, binned by two wavelength elements) 
is noisier yet faithfully tracks the SED of the low- res spec- 
trum. According to SMART, the low spectra achieve a 
signal-to-noise ratio ~ 100-300 from 9 through 15 iim 
and above 5-10 a at all wavelengths (bottom panel). 
The signal-to-noise ratio for the high-res spectrum ranges 
from ^ 50 at 10 fim to ^ 10 at 19 /im. 

To provide a check on the signal-to-noise and relia- 
bility of our low-res spectrum, we performed two tests. 
First, we extracted spectra from the two nod positions 
separately and then compared the flux densities at each 
wavelength (Figure [2]). Through 15 the spectra show 
very good agreement. The spectra become noisier at 15- 
30 /im and then produce highly uncertain fluxes at ~ 
30-38 /iTO. These comparisons allow us to empirically 
estimate the signal-to-noise per nod position as a func- 
tion of wavelength (right panel). Through 15 fim the 
median signal-to-noise is 57 per nod. From 15 /im to 
30 /im, the median signal-to-noise per nod is ~ 37. The 
median signal to noise at 30-38 /im is 6.4. The median 
signal-to-noise of the combined and extracted spectrum 
should then be ~ 81 through 15 fim, ~ 52 from 15-30 
/im, and ~ 9 at 30-38 /im. While these estimates are 
about a factor of two lower than those for the combined 
spectrum as derived by SMART, they confirm that the 
spectrum has a high signal-to-noise through ~ 30 fim. 

Second, we compared our SMART spectral extraction 
to that from the Formation and Ev olution of^lanetary 
Systems (FEPS) pipeline (jBouwman et al.ir2008[) . Fig- 
ure [3] displays the two extractions. The FEPS pipeline- 
produced spectrum consistently has a higher flux than 
the SMART extraction (top panel). However, rescal- 
ing the FEPS spectrum by 0.95 yields an almost perfect 
agreement between the two extractions (bottom panel). 
Since the FEPS and SMART-produced spectra are essen- 
tially scaled versions of one another, the differences likely 
arise from the way in which pointing offset refinements 
are treated, which compensate for the target not being 
precisely positione d at the cente r of the slit (J. Bouwman 
2010, pvt. comm., ISwain et all[20ba) . We conclude that 
our spectral extraction method is robust and interpret- 
ing the presence of spectral features, which depend on 
the relative fiux at different wavelengths, is not hindered 
by uncertainties/errors in the absolute fiux calibration 
from our extraction method. 

In part, these tests were performed to interpret strong 
features near ~ 10 fim. The IRS SLl order has a perma- 
nent bad pixel lying very close to the spectrum at one of 
the nod positions near 10 fim: failure to properly iden- 
tify and mitigate this pixel can lead to a spurious strong 
emission peak (G. Bryden 2009, pvt. comm.). Given the 
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excellent agreement between spectra extracted from dif- 
ferent nod positions and that extracted using the FEPS 
pipeline, we conclude that we have succesfuUy corrected 
for any such bad pixel-induced spectral artifacts. 

3. ANALYSIS 

3.1. Basic Properties of the EF Cha Disk: Disk 
Evolutionary State and Dust Temperature 

First, we examine the spectra, looking for any evidence 
of nebular gas emission expected if the EF Cha disk is a 
"remnant" protoplanetary disk instead of a debris disk. 
In neither the low nor the high-resolution spectra iden- 
tify evidence for H2 gas emission at 12.1 /Ltm, 17 yitm, or 
28 u rn expected for disks with c o pious amounts of gas 
(e.g. iGorti and HoUcnbach 2004; Pascu cci efal] 120061: 
ICarmona et a l. ,20Qa; Martin- Z a idi et al..l2010D . Neither 
do we find evidence for [Ne II] emission at ^ 12.8 /im, 
indi cative of trace amoun ts (~ 10~^ Mj) of circumstellar 
gas (iPascucci et al.|[2007l ). There may be a very faint, low 
signal-to-noise emission feature at ^ 6 correspond- 
ing to water gas, though its significance is less than 3 a. 
Since this emission, if real, can be produced by evapo- 
rating icy bodies it does not imply the presence of sub- 
stantial primordial nebular gas. 

Detecting gas emission lines like H2 and [Ne II] 
has proved dif ficult for younger st ars with masses like 
EF Cha's fe.g. lJuhasz et al.ll2010D . some of which like 
AB Aur have dir ect evidence for circumstellar gas 
(jBitner et al.ll2007l ). Because clear mid-IR detections of 
gas may only be possible for small fractions of the Her- 
big AeBe disk lifetime (e.g. lMartin-Zaidi et alll2010[ ). the 
absence of these lines does not guarantee the absence of 
gas in the system. However, less than ^ 1-2 % of stars as 
old as EF Cha show evidence for gas (e.g. ICurrie et al.l 
l2007d ). and most protoplanetary d isks around M ft! 2 
M(D s tars disappear prior to 3 Myr (ICurrie and KenvonI 
|2009| ) . T he EF Cha disk is also optically thin at all wave- 
lengths (jRhee et al.l 120071 ) and its emission can be well 
fit with one or two single temper ature blackblodies con- 
sistent with narrow debris belts ()Rhee et al.ll2007l Sec- 
tions 3.2 and 3.3 below). Protoplanetary disks, on the 
other hand, are typically optically thick in the mid-IR 
and are better fit by disks with large spatial extents. 
Both the low resolution and high-resolution spectra con- 
tain a strong, broad 10 fim silicate emission feature that 
is due to the presence of many small, micron-sized dust 
grains. Thus, to first order the mid-IR spectrum of EF 
Cha is consistent with tha t expected f o r a ga s-poor, dusty 
debris disk, as arg ued bv lRhee et all (|2007t) . 

Next, we construct a full SED of the EF Cha system to 
determine the temperature of the dust producing the ob- 
served IRS excess emission. To our Spitzer data, we add 
B and V data from the Tycho-II c atalog and JIIK^ data 
from the 2MASS All-Sky Survey (iSkrutskie et al.l [20061 
Figure H]) . We fit the optical data to the appropriate 
Kurucz stellar atmosphere model f or an A9 star, usin g 
the effective temperature scale from ICurrie "eTalllMl . 
Removing the stellar photosphere model from the total 
IRS spectrum, we find disk excess emission rising above 
the noise level at ~ 5-6 /^to and dominating the system 
emission by 10 //to. 

To make an initial determination of the typical dust 
temperature, we fit the observed SED to a simple two- 
temperature blackbody model, with dust temperature 



ranges of 50-400 K for the cool component and 300-1000 
K for the hot component. In doing this, we required 
that the model reproduce the observed SED in regions 
lacking obvious solid state features (A ~ 9-12 /im, 15- 
25 /im). The SED is best modeled by dust belts with 
temperatures of ^ 325 K and 490 K, corresponding to 
dust in the terrestrial zone of the system. More gener- 
ally, the SED at 30-35 ^im is consistent with a Fi, oc 
decline expected for warm dust in the Rayleigh-Jeans 
limit. Assuming a luminosity of 9.7 Lq for EF Cha and 
the blackbody relations between dust temperature and 
location - Tdust ~ 280 Kx(r/1 AlJ)-^/^{l.star /l^Sunf'^ 
these temperatures imply the presence of dust belts at 
1.0 AU and 2.4 AU. As we were able to fit the entirety 
of the observed fiux to within 20%, we conclude, based 
on these simple fits, that EF Cha lacks evidence for any 
cold debris component. 

Since the EF Cha disk spectrum does not resemble 
a featureles s blackbody as in most debris disks (e.g. 
IChen et al. ''2006). accurately determining dust temper- 
atures is complicated by the presence of solid state fea- 
tures. Moreover, accounting for these features may shift 
the best-fit dust temperatures and thus the inferred lo- 
cation of the dust. Nevertheless, we can conclude that 
EF Cha's disk emission is consistent with emerging from 
one or two debris belts at temperatures higher than the 
torr water-ice condensation temperature (~ 170 K; ~ 
8.5 AU from the star at LTE). Thus, EF Cha shows 
evidence for a warm debris disk with emission extend- 
ing from the terrestrial zone to regions analogous to the 
solar system's ast eroid belt, confirming earlier work by 
iRhee et al.1 (|200l . 

3.2. Comparisons with Protoplanetary Disk Spectra and 
Typical Debris Disk Spectra: Empirical Limits on 
Dust Composition 

The mid-IR spectrum of the EF Cha debris disk shows 
evidence for many solid state features. To investigate the 
composition of the dust disk, we first compare its emis- 
sion to that from protoplanetary disks around Herbig 
AeBe stars and typical debris disks around 1-3 Mq stars. 
Herbig AeBe star spectra exhibit a wide range of mid-IR 
spectral features, ranging from the comet-like spectrum 
of HD 1 00546 to the 10 um silicate-free spect rum of HD 
169142 (|Meeus et al.ll200ll : iGradv et al.ll2005h . 

The EF Cha spectrum is significantly different than 
that of the Herbig AeBe star HD 100546 (Figure [5]) pri- 
marily in that it lacks strong PAH emission features at 
6.3 /iTO, 7.7 /im, and 8.6 //to. Cor npared to other pro - 
toplanetary disk spectra studied by iMeeus eTall (l200l . 
EF Cha's disk shares the greatest similarity with disks 
hke HD 150193 and HD 163296 labeled as "Group Ila", 
which have broad, "boxy" solid state features at ~ 10 
microns and generally lack clear PAH emission features 
at 6.2, 7.7, 8.6, and 11.2 ^m. EF Cha shows a secondary 
bump at 11.2 /ito - comparable in wavelength to a PAH 
emission feature such as that exhibited by HD 100546 
but also easily explained by emission from olivine grains 
- lacks PAH emission features at shorter wavelengths. 

Most debris disk spectra are featureless; for example, 
even when the dust is warm enough to emit significantly 
at 10 /fm, debris disks usuall y lack the intrinsically strong 
10 /iTO silicate featu re (e.g. iDahm and Carpenter! 120091 : 
iMorales et al.l l2009l ) . The lack of silicate features indi- 
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cates tha t the dust grains a re large with diameters > 10 
(e.g. IChen et all 120061 ). In sharp contrast, the EF 
Cha spectrum shows a complex 10 fim peak and spectral 
structure longwards of « 15 fim, indicating that the disk 
contains many grains smaller than 10 fim. This behav- 
ior is qualitatively similar to that of other warm, high- 
luminosity debris disks around young stars, s uch as those 
found around HD 11 3766 A and HD 172555 (|Lisse et al.l 
[TOOS, 2009, Figure [5]) . 

High mass fractions of crystalline pyroxene can also 
be ruled out since they have multiple, narrow emission 
peaks in the 8-13 range resulting in a spectrum far 
less smooth than observed. Likewise, the spectrum can- 
not contain as high of an abundance of crystalline olivine 
as that exhibited by comets and some Herbigs like the 
HD 100546 disk: such a high abundance would result in a 
narrow maximu m at ^ 9.7 iim and a higher maximum at 
^11.2 fim (see iLisse et al]|2007al ). the opposite of what 
is observed. The general shape of the EF Cha spectrum is 
instead consistent with a high abundance of amorphous 
silicates and a secondary population of olivine. Amor- 
phous silicates produce much broader 8-13 fim peaks, 
resulting in a secondary maxima at ~ 17-19 fim, and 
are abundant in Group Ha objects (e.g. iBouwman et al.l 
1200 It iLisse et al.ll2007aD . However, superimposed upon 
this shape are peaks at 10, 11, 19, 23, 28, and 33 fim. 
The existence of these peaks indicates that additional 
species besides amorphous silicates are present. 

3.3. Comparisons with Mineralogical Models 
3.3.1. Description of Mineralogical Model 

As shown by the empirical comparisons, the compo- 
sition of mid-IR emitting dust from EF Cha's debris 
disk is dominated by amorphous silicates but must in- 
clude high mass fractions of other species. To identify 
these additional species, we fit the EF Cha IRS spectrum 
with a detailed mineralogical model based on laboratory 
thermal inf r ared emission s pectra described in detail by 
ILisse et al.l (l200l [2007allbl . l2008ll . Below, we reiterate 
the key components of the model and describe its appli- 
cation to our data. 

Thermal Dust Emission - The emission from dis- 
tribution of dust is given by 



emission spectra drawn from measurements for ran- 
domly oriented, /im-sized powders, which yield Qabs 
for each species. Sources for the emission spec- 
tra are d etailed in t h e Sup plemental Online Mate- 
rial from ILisse et aLl (|2006[ ) and among others in- 
clude contributions from the Jena spectral library 
( |http://www.astro . u ni-iena.de/L a borato r y/OC DB) and 
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where T is the particle temperature for a particle of ra- 
dius a and composition i at a distance r from the cen- 
tral, the distance from Spitzer to the dust is A. B>, is the 
blackbody radiance at wavelength A, Qabs is the emission 
efficiency of the particle of composition i at wavelength 
A, dn/da is the differential particle size distribution of 
the dust. The sum is over all species and all particle 
sizes. Our spectral analysis consists of calculating the 
fiux for a model distribution of dust and comparing the 
result to the observed fiux. The predicted flux from the 
model depends on the 1) dust composition (affects the 
location of spectral features), 2) particle sizes (affects 
feature-to-continuum contrast), and 3) the particle tem- 
perature (affects strength of short versus long wavelength 
features). 

Dust Composition- We use thermal laboratory 



from iKoiki (l2000l) : iKoike et al.l (l2002l ): Ichihara et aTl 
()2002| ) for silicates. The library of material spectra were 
selected by their reported presence in interplanetary dust 
particles, meteorites, in situ comet me asurements , proto- 
planetary disks, and debris disks f see ILisse et al.l [20061. 

In total, the library includes spectra from 80 different 
species. The list of materials considered include multi- 
ple silicates in the olivine and pyroxene class (forsterite, 
fayalite, clino- and ortho-enstatite, augite, anorthite, 
bronzite, diopside, and ferrosilite) . The model also 
includes phyllosilicates (such as saponite, serpentine, 
smectite, montmorillonite, and chlorite); sulfates (such 
as gypsum, ferrosulfate, and magnesium sulfate); ox- 
ides (including various aluminas, spinels, hibonite, mag- 
netite, and hematite) Mg/Fe sulfldes (including pyrro- 
htite, troilite, pyrite, and niningerite) ; carbonate min- 
erals (including calcite, aragonite, dolomite, magnesite, 
and siderite); water-ice, clean and with carbon dioxide, 
carbon monoxide, methane, and ammonia clathrates; 
carbon dioxide ice; graphitic and amorphous carbon; and 
neutral and ionized polyaromatic hydrocarbons (PAHs). 

Dust Size Effects - We consider a dust size range 
of 0.1 to 1000 fim in the model flts, with sharp emission 
features arising mainly from particles of 0.1 - 10 fim, and 
baseline continuum mainly from 10 - 1000 fim particles. 
The particle size effects the dust emissivity as 

1 - Emissivity{a, A) = [1 - Emissivity{lfim, X^°-/^t^"^ 

(2) 

The particle size distribution (PSD) is fit at log steps in 
radius, i.e., at 0.1, 0.2, 0.5, 1, 2, 5, ...10, 20, 50 jim, etc. 
Particles of the smallest sizes have emission spectra with 
very sharp features, and little continuum emission. Par- 
ticles of larger sizes - greater than 10 fim - tend to 
be optically thick, and contribute only featureless con- 
tinuum emission. Particles larger than 1000 fim are not 
directly identifiable from our modeling and contribute no 
emission. We assume the same particle size distribution 
daioY each species. 

Dust Temperature - Dust particle temperatures are 
determined at the same log steps in grain radius used to 
determine the PSD. Dust temperature is a function of 
dust composition, size, and astrocentric distance. The 
highest temperature for the smallest particle size of each 
species is varied freely and is determined by the best fit to 
the data. The largest, optically thick particles (1000 fim) 
are set to the LTE temperature, and the temperature 
of particles of intermediate sizes is interpolated between 
these two extremes by radiative energy balance. 

Accuracy /Reliability of the Model and Fitting 
Method - This model has been successfully applied 
to both solar system objects and protoplanetary/debris 
disks around other stars. In particular, the model ac- 
counts for the emission spectra of come t 9P/Tempel 1 
observed during the De ep Impact mission (lA'Hearn et al.l 
I2005t ILisse et al.ll2006[) . Prior to the collision with the 
spacecraft, the comet's IR spectrum was near feature- 
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less, consistent with a coma dominated by large black- 
body grains, similar to debris disk spectra around most 
stars. After collision, the spectrum showed more than 16 
pronounced emission features, strong evidence for » 10'' 
kg of small, submicron to micron-sized grains consisting 
of amorphous and crystalline silicates, amorphous car- 
bon, carbonates, phyllosilicates, PAHs, water gas and 
ice, and sulfides (jLisse et al.ll2Q06[ ). As evidenced by the 
copious presence of fragile, volatile water ice grains in 
the ejecta, the differences in the pre- and post-impact 
spectrum must be due to simple de-aggregation of weak, 
large 100 /im to 1 cm), optically thick fractal dust 
particles in the delivered shock wave, rather than miner- 
alogical transformation of the dust at high temperatures 
and pressures due to the delivered energy. 

Direct compositional measurements from cometary 
sources also provide strong evidence that our modeling 
approach 1) is accurate for dust in the solar system and 
2) can infer the composition of dust around other stars. 
A comparison of the derived Tempel 1 material com- 
position with the materials found in the STARDUST 
sample return show good agreement (Hanner and Zolen- 
sky 2010). Moreover, the direct measurements of species 
within interplanetary dust particles (IDPs) of cometary 
origin, in particular the GEMS (glass with embedded 
metal and sulfides) grains located within the matrices 
of "chondritic porous" IDPs, provide good evidence that 
cometary/meteoritic provide a good spectral match to 
"astronomical" silicates such as those infer red from spec- 
tra of protoplanetary and debris disks ([Bradlev et al.l 
[1999). 

For both Tempel 1 and debris disks around other stars 
like EF Cha, the dust responsible for mid-IR solid state 
features is produced from collisions and must be freshly 
produced as it is small and easily removed from the sys- 
tem by radiation pressure, Poynting- Robertson drag, and 
stellar wind drag. Thus, applying our model to debris 
spectra from other stars can 1) identify the chemical com- 
position of parent bodies from which the debris originate 
and 2) determine the sizes of debris particles, which may 
identify whether the debris must originate from large 
scale, catastrophic ii npacts and 3 ) determine the range of 
dust temperatures. iLisse et all (|2008l I2OO90 argue that 
debris emission from disks surrounding HD 113766A and 
HD 172555 originates from debris created by large-scale 
coUisional grinding/giant impacts associated with terres- 
trial planet formation and/or stochastic events such as 
that responsible for forming the Moon. Likewise, we here 
apply the model to the EF Cha spectrum to constrain the 
composition of parent bodies and investigate the origin 
of the debris. 

Our method differs from other widely-used methods in 
that it uses physically plausible emission measures from 
randomly-oriented powders rather than theoretically- 
derived Mie values. The free parameters of the model are 
the relative abundance of each detected mineral species, 
the temperature of the smallest dust particle of each 
species, and the value of the particle size distribution. 
Best-fits are found by a direct search through phase space 
- composition, temperature, and size distribution. 

3.3.2. Results of Fits From Our Mineralogical Model 

Figure [S] shows our fit to the EF Cha flux, presented in 
emissivity space (spectrum divided by an average Planck 



function of the entire dust population) for clarity, with 
the emissivity of individual constituent species labeled 
and shown as colored lines. Table 1 lists the relative 
abundances of each species. The best-fit model requires 
the presence of two warm dust belts, one with micron- 
sized dust at 600 K (LTE at - 430K, or - 1.3 AU) 
and one with large dust grains (> 20 ^m) at an LTE 
temperature of ~ 220 K (~ 5 AU). 

The best-fit model successfully reproduces the entire 
IRS spectrum and yields a small per degree of free- 
dom (x^/dof — 1.06). To provide a statistical measure 
of the goodness-of-fit, we determine the 95% confidence 
limit given the number of degrees of fr eedom for a given 
reduced fsee lLisse et"alll2008l [2009t ). The total num- 
ber of free parameters - including the number of rel- 
ative abundances, number of hottest particle tempera- 
tures, and power law index - is 32. Given the number of 
spectral points fit (272), the number of degrees of free- 
dom is 240. Based on these values, there is a 95% chance 
a model with reduced less than 1.15 is a good pre- 
dictor of the Spitzer data. Given the excellent fit of the 
model from 5 to 38 /im, it is plausible that the chemical 
abundances included in the model comprise the domi- 
nant detectable species in the disk. 

As suggested by the empirical comparisons, the EF 
Cha spectrum is olivine rich, containing large relative 
molar abundances of amorphous olivine and forsterite 
i^moies.rei = 0.73 and 0.46, respectively) responsible for 
the large width of the 10 urn feature and the presence of 
a secondary peak at ~ 11 /im. Fayalite also is present in 
low abundance. The model x^i^ is 1.07 with this species 
removed, lower than the 95% confidence limit, so evi- 
dence for this species is very weak. The emissivity slope 
and features at ^ 20-35 /xm indicate the presence of wa- 
ter ice and metal sulfides, which are also abundant in the 
EF Cha disk {J^moies,rei — 0.44 and 1.45, respectively)). 
The spectrum contains a small abundance of pyroxene 
but lacks any clear evidence for water gas, carbonate, or 
PAH emission typical of primitive cometary material. 

Interestingly, the EF Cha debris disk appears to con- 
tain a very high abundance of phyllosilicate species - 
smectite, talc, and saponite - responsible for the ampli- 
tude and shape of the pronounced peak at 10 /xm, 15 
and 20-23 {^moies,rei — 0.2 total). Figure[71 display- 
ing disk emissivity with silicate contributions subtracted 
out, more clearly shows the importance of phyllosilicates 
and other major constituents like metal sulfides. In par- 
ticular, phyllosilicates have singular sharp features and 
are required to reproduce the strong residual 10 /xm emis- 
sion and complicated slope of the emissivity at 17-25 ^m. 
SiO and silica have peaks at 8.8 to 9.2 /xm wavelengths, 
much shorter than 10 /im, and their inclusion in models 
clearly does not improve the x^ parameter. Without the 
presence of phyllosihcate species, the x^j^ of the spectral 
fit is far worse - 3.1 for the removal of smectite alone. 
The range of abundances from any model passing the 
95% confidence limit is narrow, typically ~ 10% from 
the best-fit value. 

Given the numerous strong solid state features, the 
EF Cha disk must contain large numbers of small grains. 
The best-fit model includes a minimum grain size of 0.1 
/im: two orders of magnitude smaller than grain sizes 
for most IRS-studied debris disks. Moreover, the par- 
ticle size distribution for the best-fit model is dn/da oc 
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g^-4.o_ 'pj^ig steep distribution implies the presence of 
many more fine particles than for a collisional equilib- 
rium size distribution (dn/da oc a~^'^). 

While we can only directly sense the presence of grains 
with sizes < 20 /im, grains up t o millimeters in size have 
been found in debris disks dManess et al.l 120081 ) . More- 
over, micron-sized dust results from the collisional cas- 
cade induced from colliding, km-sized planetesimals and 
thus implies the presence o f much larger bod i es (e.g. 
Kenvon and BromlevI 120081: IWvatt et all 120081: iKrivovl 
20101 ). The size distribution of these fragments are 
predicted from theory to follow a power law behavior 
roughly continuous with that for directly detec table de- 
bris dust (e.g. see iKenvon and BromlevI [20081 ) . There- 
fore, we can get a rough estimate of the total mass of 
fragments involved in the collisions by extending the par- 
ticle size distribution up to larger sizes of order of 1 mm. 

Assuming a maximum grain size of 1000 /ito, we esti- 
mate a total dust mass of ^ 1.5x10^'^ g, comparable to 
the mass of an asteroid with a bulk density of 2.5 g cm~'^ 
and radius of ~ 110 km. This parent size is a lower limit 
since it is based on the total mass of grains with sizes less 
than 1 mm. Larger, 100-1000 km-sized bodies are likely 
present since they are required to initiate and sustain 
the collisional cascade that produces the < 1 mm-sized 
debris (e.g. iKenyon and Bromley 2008, 2010). 

The large number of small grains must require a replen- 
ishment source since they would o therwise be rapidly re- 
moved by radiation pre ssure (e.g. iBackman and Paresca 
119931 : IChen et al.ll2006D . Assuming a stellar mass of 1.9 
Mq, a stellar luminosity of 9.7 L©, and a grain density of 
3 g cm~^, the blowout size l imit for spherical gra i ns in EF 
Cha, using the equation in lCurrie and KenvonI ()2009l ) is 
~ 2 /im. Thus, the EF Cha debris disk contains a high 
mass of dust dominated by particles comparable to or 
smaller than those dominating typical debris disks and 
also far smaller than the nominal blowout size limit for 
the system. Along with the stabilty of the 1983 - 2007 
IRAS/MSX/Spitzer photometry, the copious amount of 
small dust indicates that there must be a source for small 
particle replenishment. 

3.3.3. Tests of Phyllo silicate Identification 

Additional tests are required to more reliably constrain 
the presence of individual species responsible for the 
EF Cha debris emission, especially phyllosilicate species. 
First, although our best-fit model includes phyllosilicate 
species our solution is not unique if a phyllosilicate-free 
model reproduces the observed spectrum within the con- 
fidence limit we adopt to identify good-fitting models 
(the 95% confidence limit). Second, the signatures from 
phyllosilicate species and other species like metal sulfides 
are weaker than the amorphous silicates also present in 
the disk. Since they are weaker their presence is more 
likely to be falsely identified in real spectra with random 
noise patterns. 

To address these issues, we carried out two additional 
sets of model fits. First, we reran the fits with all phyl- 
losilicate species removed, in order to determine whether 
all good-fitting models require their presence. Second, we 
fit our models to two independent EF Cha spectra ob- 
tained concurrently at separate IRS nod positions. Com- 
paring the two model results constrains the sensitivity 
of our identification of different species to photon noise. 



Moreover, since the spectra at different nod positions are 
extracted from different regions on the IRS detector, fits 
to separate nod positions test the sensitivity of our identi- 
fications to different IRS systematic noise contributions. 

Figure m displays the best- fit model obtained with con- 
tributions from all phyllosilicate species removed. If the 
best-fit model fares only slightly worse in reproducing 
the inferred emissivity with and without silicate contri- 
butions, then the EF Cha disk composition is underdeter- 
mined and we cannot constrain the existence and mass 
fractions of phyllosilicate species. As shown by Figure 
m the of the best-fit model lacking phyllosilicates is 
substantially worse than for the model including phyl- 
losilicates. While the general shape of the emissivity is 
accurately reproduced, the phyllosilicate-free model pre- 
dicts a strong peak at ~ 11.25 /im, which is not observed, 
and incorrectly predicts the flux ratios of peaks at 10 fj,m 
and 11 fj,m. The model also significantly underpredicts 
the flux at 23 /im. 

The silicate-subtracted emissivity even more strikingly 
reveals the failure of our phyllosilicate-free model (Fig- 
ure [9]). The model emissivity is relatively flat from 8 fj,m 
to 11 ^m. With the silicate contribution from the model 
subtracted out, the observed emissivity has local maxima 
at ~ 9.5-10 fim and 23 fim well above that of the model. 
Local minima at 8-9 fim and 11-12 fim lie well below 
that of the model. The model also lacks the structure 
apparent at ^ 14-20 fim in the silicate-subtracted emis- 
sivity from the data. Compared to the phyllosilicate- 
free model, the best-fit model including phyllosilicates 
clearly reproduces the emissivity far better. Quantita- 
tively, the phyllosilicate-free model can be ruled out as 
a good-fitting model. It has a reduced value ~ 40% 
larger (x^/dof = 1.47 for 240 degrees of freedom). This 
difference in is statistically significant, since the prob- 
ability that the phyllosilicate-free model fits the Spitzer 
data is ~ 0.0002. 

Table [5] summarizes the modeling results for our fits 
to the separate nod positions. The number of moles and 
weighted surface area for the species as a whole show 
exceptional agreement with the results derived for the 
combined spectrum. While a few of the abundances do 
vary - e.g. for amorphous olivine, amorphous pyroxene, 
smectite - they do so by less than the 10-20% systematic 
uncertainties on the abundances in all cases. While only 
one nod position requires the presence of a small amount 
of orthoenstatite, the presence/absence of other species 
is completely consistent from one nod position to the 
other. 

Most importantly, fits to both nod positions require 
the presence of the three phyllosilicate species identified 
from fitting the combined spectrum. The total weighted 
surface area for these species varies little between nod 
positions (0.32 vs. 0.34), and as a result the molar abun- 
dances vary little (0.19 vs. 0.20). Thus, our quantitative 
estimate of the phyllosilicate contribution to the EF Cha 
debris disk does not depend on the nod position used. 
Since fitting to different nod positions probes the sensi- 
tivity of our results to noise characteristics, our results 
show that our identification and abundance estimate of 
phyllosilicate species is robust against the infiuence of 
random and systematic noise. 

At a minimum, we can conclude that some suite of 
species whose emissivities are extremely similar to the 



7 



laboratory-measured emissivities for our suite of finely- 
powdered phyllosilicate species are abundantly present 
in the EF Cha debris disk. Given the complexity of the 
phyllosilicate family, and the large range of variants for a 
given phyllosilicate (e.g., the intechangability of Ca and 
Al, as well as Fe and Mg) it is not totally surprising that 
we may not yet have the exact laboratory match for the 
EF Cha spectral signature. Moreover, the precise shape 
of the spectral features depend on the exact details of 
how the minerals in the parent bodies are fragmented, 
combined, heated and cooled. A priori, we cannot rule 
out the existence of species not included in our chemical 
model that can mimick the emissivity of phyllosilicate 
species. These uncertainties preclude any proof that the 
EF Cha disk must contain phyllosilicate material. 

However, the presence of abundant phyllosilicate ma- 
terial in a variety of solar system contexts - e.g. car- 
bonaceous chondritic meteorites, Comet 9P/Tempel- 
1, main belt asteroids jT omcoka and Ohnishi 2010; 
Morlok ct al. 2010; Lissc et_aL| i2006t iVilas and GaffevI 
1989; Bus and Binzcl 2002)- lends powerful support to 
our inference that the species in question is indeed phyl- 
losilicate species. If the species in question is not phyl- 
losilicates, chemical processes must have been able to 
remove these species selectively and thoroughly, which 
seems implausible given the wealth of other species 
present. Given the extensive library of species used 
for our spectral f itting and rejected as poorly matching 
(|Lisse et al.ll20Q6[ ). it also seems unlikely that some other 
species is mimicking the phyllosilicate species' signatures. 

3.4. Comparisons to the Mineralogy, Dust Sizes, and 
Atomic Abundances of Comet Tempel-1 Ejecta and 
Warm, High- Luminosity Debris Disks 

To provide an additional framework for interpreting 
our results, here we compare the EF Cha dust mineral- 
ogy, dust sizes and atomic abundances to those for the 
Comet Tempel-1 ejecta and warm, young debris disks we 
have m odeled in previou s papers: HP 1137 66A and HD 
172555 (" Lisse et al.l[2006l [20Q7aI [2Q0l [2009h all of which 
demonstrate strong mid-IR features like EF Cha. For 
both Comet Tempel-1 and HD 113766A, the 9-12 /xm 
flux exhibits a broad plateau with two peaks of equal 
amplitude at 10 /xm and 11 /im due to abundant 
forsterite, enstatite, ferrosilite, and (to a lesser extent) 
amorphous olivine and pyroxene. Best-fit models for the 
Comet Tempel-1 ejecta and HD 113766A also require 
the phyllosilicate species smectite nontronite. While HD 
172555. has a single sharp peak, it is at ~ 9 iim with 
a gradual decline in flux from 10 to 12 fim dLisse et al.l 
[2009) . Our modeling identifies this peak as due to a com- 
bination of amorphous silica, tektite, and SiO gas. EF 
Cha's spectrum is different from all of these, exhibiting 
a strong 10 jim phyllosilicate peak superimposed upon a 
smoother 10 fim amorphous silicate feature and the 9.3, 
9.8, and 11.2 fj,m sharp features of crystalline olivine & 
pyroxene. The phyllosilicate molar abundance for non- 
tronite, talc, and saponite combined is far higher than 
that for HD 113766A or Comet Tempel-1 (Nmoies,rei = 
0.20 vs. 0.03 and 0.07) and thus much more significant. 

The abundant small grain population for EF Cha is 
most consistent with a large mass fraction of debris 
emission being actively produced, either by a single, 
massive catastrophic collision followed by grinding of 



the fragments or an abnormally high, steady collision 
rate. The agreement in flux density between the 1983 
IRAS, 1996 MSX, and 2007 Spitzer data is at least con- 
sistent with the emission being maintained at a high 
level for 20 years or more. Timescales for grain re- 
moval by Poynting-Robertson drag are longer than this 
limit fe.g.lBackman and Parescel ll99l lChen et al.l l2006l: 
ICurrie and KenvonI l200'9l )~ however, removal timescales 
for the smallest grains due to radiation pressure are more 
comparable (e.g. orbital timescales of 1-100 yr), indicat- 
ing that at least the smallest grains may require active 
replenishment. Similarly, the presence of small grains 
and very steep particle size distribution is qualitatively 
consistent with that seen in the Tempel 1 ejecta produced 
by the Deep Impact collision and with fine dust created 
by a hypervelocity impact found around HD 172555. 

Like HD 11 3766 A and Comet Tempel-1, the EF Cha 
spectrum may indicate the presence of water ice. HD 
113766A has two debris belts beyond the ice line and a 
warm one inside of it. Our spectral modeling suggests 
that EF Cha contains a warm belt well interior to the 
ice line and second one close to the ice line, at 220 K at 
7.6 AU. Water ice present in the interior EF Cha belt, in 
a radiation field with a local equilibrium temperature of 

600 K, should be very quickly vaporized and thus not 
present in the mid-IR spectrum unless it has been freshly 
ejected from a parent body. Collisions producing the 
observed water ice then either must have occurred very 
recently or must be occurring in the outer, 220 K belt. 
The belt at 7.6 AU may even provide a reservoir for some 
icy dust concentrated in the warm, terrestrial zone belt. 
In either case, the presence of water ice requires an active 
replenishment source. Moreover, if water ice is being 
ejected from parent bodies in the 7.6 AU region near 
the ice line, the EF Cha spectrum may show evidence for 
faint, water gas emission. If water ice is indeed present 
there must be an active replenishment and a continual 
release of gaseous water. We do not detect water gas at 
a statistically significant level (see Section 3.1), although 
the expected amount is difficult to quantify. Thus, we do 
not rule out its presence at the observed signal-to-noise 
level. 

Within the context established by data for comets, me- 
teorites, protoplanetary disks, and warm debris disks, 
EF Cha's olivine mass fraction indicates that the sys- 
tem bridges a gap between relatively primitive unpro- 
cessed solids and heavily processed bodies (Figure[Tni top 
panel). Primitive bodies coming in contact with water, 
and thus aqueously altered, gain in their mass fractions 
of olivine and olivine-to-pyroxene ratios^. Differentiated 
bodies and carbon-rich bodies processed at high temper- 
atures have an even higher abundance of olivine and lie in 
the upper- right hand corner of the plot. Samples with a 
low percentage and relative fraction of olivine occupy the 
lower- left corner and include comets (SW3, Holmes, Wild 
2, and Tempel 1), and prot oplanetary disks ( HD 100546, 
HD 9 8800, HD 163296) (iSitko et all \200§. iLisse et al.l 
[2001. [2007al [2008.1 1. Those in the upper-right are con- 
sistent with differentiated and/or pyrolised bodies in- 

^ Here we include phyllosilicates in the pyroxene mass frac- 
tion, assuming these species a re formed by alteration of pyroxene 
by warm, rea ctive wa ter (s ee ITomeokal 119901 : IMorlok et al.l 120101 : 
INagahara and Ozawail2011l') . 
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elude ancient debris disks surrounding white dwarfs (GD- 
362 G29-38) and the centaur SW-1 (IStansberry et aLl 
[200I iReach et aLll2009bl: [Jural [200I lJura et all 12009^ 

In this regard, EF Cha is roughly halfway between Hale- 
Bopp/Tempel-l and HD 113766, which respectively iden- 
tify the boundaries of the the primitive, unprocessed 
group and the highly processed, differentiated/pyrolysed 
group. 

The bottom panel of Figure [TUl compares atomic abun- 
dances for the EF Cha debris disk to that for the Sun 
and for other debr is disks, p r otopl anetary disks and 
comets (e.g. see iLisse et all l2008| ). Like other de- 
bris/protoplanetary disks and comets - HD 172555, 77 
Corv, HD 100546, etc - the EF Cha emitting debris 
has heavily subsolar abundances of hydrogen, carbon and 
oxygen. EF Cha has solar abundances of the main refrac- 
tory elements Si, Mg, Fe, and Al, a subsolar abundance of 
calcium, and supersolar abundances of sulphur. Overall, 
EF Cha's atomic abundance pattern shows the strongest 
similarity with that of HD 113766A, another luminous 
warm debris disk of comparable age whose spectrum 
shows evidence for massive terrestrial-zone collisions. To 
first order, EF Cha's debris dust may be thought of as 
a more unprocessed/amorphous, water-rich analogue to 
HD 113766A's. 

3.5. Interpretation of Modeling Results 

Phyllosilicate species are tracers for the past pres- 
ence of warm reactive water, since they are pr oduced by 
the aqueous alteration of anhydrous rock (e.g. iTomeokal 
[1990). Based on our modeling results, the emitting dust 
from EF Cha's debris disk plausibly includes a high abun- 
dance of phyllosilicate grains and smaller abundances of 
water ice and metal sulfides in the terrestrial zone re- 
leased from parent bodies by collisions. Because phyl- 
losilicates shoul d slowly decompose at temperatures o f 
T - 600-700 K (|Akailll992t[Nagahara and Ozawal[20H . 
they can exist within much of the terrestrial zone on the 
order of years after release from their parent body. How- 
ever, stable low pressure phyllosilicate formation over 
millions of years, whether from solar nebula condensa- 
tion or from aqueous alteration in the interior of an aster- 
oidal water-rich body, requires bodies that were formed 
at temperatures < 300 K (e.g. iFeglev and PrinnI 119891 : 
iNagahara and Ozawall201l[ ). Just as the presence of icy 
parent bodies scattered into the terrestrial zone from be- 
yond the ice line that can be inferred from water ice 
spectral features, we can conclude that the current warm 
dust belt arises in part from solids imported from cold 
(< 300K) regions, more distant from the star. The com- 
position of the parent bodies responsible for EF Cha's 
debris emission is then likely not characteristic of plan- 
etesimals that, ^ 5-10 Myr prior, grew from micron to 
millimeter-sized grains in the warm, terrestrial regions 
of its protoplanetary disk. Rather, the current debris 
disk composition arises in part from solids imported from 
colder regions, more distant from the primary star. 

The large abundance of small grains, evidence for high 
collision rates, and atomic abundance pattern makes EF 
Cha's debris disk comparable to the high-luminosity, 
warm debris disks around HD 113766A. However, its 
chemical composition is very different, given its high 
phyllosilicate abundance and moderately low pyroxene 
abundance consistent with aqueously altered primitive 



material. If large-scale terrestrial planet formation is oc- 
curing in EF Cha like HD 113766A, then it is occuring 
from a feedstock of water-rich C- or D-type asteroids, 
not water-poor S-type asteroids as inferred by Lisse et 
al. (2008) for HD113766A. Combining our disk chem- 
istry and grain size results, we conclude that the EF Cha 
debris disk traces 1) a high level of terrestrial zone colli- 
sions between bodies which, at least in part, 2) originated 
from beyond the terrestrial zone, and 3) have processed 
phyllosilicates from primordial silicates at some time in 
the past. 

4. DISCUSSION 

4.1. The Importance of EF Cha's Phyllosilicate Species 

The detection of abundant phyllosilicates in EF Chas 
debris disk is of particular interest to studies of the early 
solar system and the habitability of terrestrial planets. 
Phyllosilicate formation results from the reaction be- 
tween anhydrous rock and water and thus its presence 
in a debris disk requires an abundant source of water- 
rich planetesimals. Meteoritic data shows evidence for 
numerous "fine-grained" phyllosilicate rims surrounding 
chondrules and calcium-aluminum-rich inclusions in car- 
bonaceous chondrites and meteorites - i.e. the par- 
ent bodies that are a likely source for Earth's water 
(e.g. [Robert et all [2000t) - with abundances similar to 
phyllosilicat e-rich asteroids in the outer asteroid belt (r 
> 2.5 AU) (|Gradie and Tedescdll98^ iBeh et al.l ll9891. 
If phyllosilicate species are detected in debris disk spec- 
tra, their presence then requires an abundant source of 
water-rich planetesimals. Thus, the presence of phyllosil- 
icates in the terrestrial zone of disks surrounding young 
stars may trace water delivery in young planetary sys- 
tems. 

Moreover, because phyllosilicates trace the existence of 
water in disks, they may be important signposts for the 
future development of primitive life. Indeed, the pres- 
ence of abundant phyllosilicates in the most ancient re- 
gions of Mars' surface has been used in part to argue 
that Mars may have once supported habitable conditions 
([Bibring et alJ[2OO60 . Surface regions rich in phyllosili- 
cates may be prime locations for finding early, primitive 
organ isms ([Farmer and Des MaraisI [19991 : [Qrofino et al.l 

mm . 

4.2. Phyllosilicate Production and Water Delivery 

Planetsimals collisions are responsible for releasing 
small grains that allow us to identify phyllosilicate 
species. However, it is important to emphasize that the 
collisions between anhydrous rock and hydrous planetes- 
imals are not likely for the production of the phyllosil- 
icate species. Impact processing tends to devolatilize 
planetesimals and should convert ai iy phyllosilicates on 
the surface back to pyroxene (e.g. iMorlok et all [201(1 
see section 4.3 of this work). Instead, phyllosilicates 
species are formed from the combination of pyroxene 
and warm, reactive water, where the water likely resulted 
from ice melted interna l heating due to ^^Al deca y (e.g. 
[Desch and Leshid [2004 [Merk and Prialnik[ [20061 , and 
references therein). Thus, it is more likely that the colli- 
sions are simply revealing aqueously altered rock as op- 
posed to being a process that actively alters anhydrous 
rock on the planetesimals' surfaces. 
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Recent studies of the solar system's main belt aster- 
oids have yielded precisely this - d etections of water and 
aqueously altered produ cts (e.g. ICampins et al.l 120101 : 
iRivkin and Emervll2010D . The total dust mass required 
to reproduce the EF Cha disk emission is comparable in 
mass to a 110 km-radius asteroid: a parent body of com- 
parable size and densit y is then massive enough to have 
contained liquid water ()Merk and Prialni3l2006f) . If our 
conclusion about the presence of phyllosilicate species is 
correct, it implies that asteroidal processing such as that 
identified from solar system asteroids may be quite ad- 
vanced by 10 Myr. 

The mechanism for the delivery of water-rich plan- 
etesimals to the terrestrial zone (needed for phyllosil- 
icate production) is less unclear, although there are 
several plausible mechanisms. After the planetesimal 
population in the outer disk is significantly depleted, 
the remaining planetesimals could be dynamically scat- 
tered inwards and incorporated into forming terres- 
trial pla nets, especially after the solar nebula has dis- 
sipated (iGoldreich et al] 120041: iFord and Chiand 120071: 
iMatsumura et al.ll2010() . Water-rich planetesimals could 
also be rapidly delivered from asteroid belt regions to 
the terrestrial zone by gas drag during the protoplane- 
tary disk phase and later incorporated into larger bodies 
(e.g. [Cicsla and Lauretta 2005)., or form nearly in situ 
if the disk is cold. In other systems, volatile-rich planets 
formed beyond the ice line could fail to become the cores 
of gas giants and migrate into the terrestr ial zone throug h 
planet-disk or planet-planet interactions (Kuchner'1 20031 ). 

From the standpoint of our mincralogical analysis, all 
of these scenarios are equivalent. EF Cha's disk is a 
debris disk, not a protoplanetary disk. Therefore, we 
do not know exactly when or how the water needed for 
phyllosilicate production was first introduced to the ter- 
restrial zone in EF Cha. 

4.3. The Detectability of Phyllosilicate Species in 
Protoplanetary Disks and Other Debris Disks 

A more answerable question is why so far only EF Cha, 
among all debris disks and protoplanetary disks, shows 
strong evidence for abundant phyllosilicates. First, de- 
tecting phyllosilicate species with Spitzer requires dust 
warm enough to effectively emit at the 10-25 /Ltm with 
sufficiently small sizes (< 10 /im) to show solid state 
features. Both models for terrestrial planet formation 
and Spitzer surveys indicate that warm debris dust 
around young stars is rare, because the processes that 
can make the dust are fleeting and/or yield too few 
small collisional fragments for detectable solid state fca- 
tures (e. g. Kcnvon and Bromlcv 2004; Chen ct al. 2006; 
ICurrie et'"al.ll2007al: fCorlova et al.ll2007t ). Thus, the sub- 
set of young debris disks that could possibly yield de- 
tectable mid-IR features of phyllosilicates is small com- 
pared to the sample of all young debris disks. 

Second, experiments on phyllosilicate-rich chondritic 
meteorites indicate that collisions cannot be too high- 
energy, else the phyllosilicates will be destroyed. 
iMorlok et al.l ()2010[ ) subjected matrix sections of the 
phyllosilicate-rich Murchison CM chondritic meteorite to 
projectile collisions with peak pressures ranging between 
10 and 49 GPa. The weakly shocked samples (-- 10 GPa) 
corresponding to low-energy collisions have spectra dom- 
inated by phyllosilicates, in particular serpentine with a 



broad 10 peak. For peak pressures between 20 and 
30 GPa, the Murchison samples retain mid-IR signatures 
of phyllosilicate species but crystalline olivine's peak at 
~ 11.2 /zm begins to emerge. At higher pressures (~ 36 
GPa) corresponding to high-energy collisions, the phyl- 
losilicate signatures are lost due to serpentine decompo- 
sition and devolatization and the spectra better resemble 
those for debris disks with weaker/absent phyllosilicate 
emission such as HD 113766A and HD 69830. Thus, 
EF Cha's debris disk plausibly represents a particularly 
rare laboratory for studying phyllosilicates in other plan- 
etary systems, where collisions are sufficient to produce 
detectable debris but not so catastrophic as to eliminate 
evidence for aqueously altered planetesimals. 

While numerous protoplanetary disks in nearby star- 
forming regions have been targeted by Spitzer, stud- 
ies typically do not incorporate phyllosilicate species 
in their spectral modeling and instead focus on amor- 
phous silicates, for s terite, enstatite, and silica (e.g. 
Sargent et all l2009l : IWatson et all 120091 : lJuhasz et all 



20101) ^. Recent model s for p rotoplanetary disk spec- 
tra bv lMorris and DeschI ()2009[ ) argue that phyllosilicates 
can be identified at ~ 10 /xm and, especially, 20-25 if 
they comprise more than 3% of the dust by surface area, 
as is the case for EF Cha, and if the dust is fine enough 
to be optically thin ^. 

However, our analysis cannot yet constrain the fre- 
quency or duration of phyllosilicate production in debris 
disks. No other debris disk has been reported to date 
as having strong, clearly phyllosilicate derived emission 
features. Further, since debris disk grains are generally 
too large to produce clear solid-state features anyway, 
we cannot conclude that phyllosilicate production in de- 
bris disks is fleeting or rare. More precisely, determin- 
ing exactly when and how often phyllosilicates form, and 
thus when water is delivered to the terrestrial zone of 
young planetary systems, will require modeling more pro- 
toplanetary disks and debris disks with solid state fea- 
tures using sophisticated mincralogical models like those 
presented here. We are encouraged that future work will 
produce more and more phyllosilicate detections, how- 
ever, given the preponderance of likely phyylosilicate rich 
C- and D-type asteroids in the outer main belt of the so- 
lar system, and the suggestion of phyllosilicates in comets 
(Lisse et al. 2006, 2007a, Reach et al. 2009a). 

Spectral modeling of protoplanetary disks is also ar- 
guably more complicated, given the effects of disk ffar- 
ing/shadowing, variability, and the inability to probe the 
disk midplane (where planetesimals should be located) if 
the disk is optically thick. However, many Herbig AeBe 
stars with protoplanetary disks in the lJuhasz et all studv 
(HD 36112, HD 37258, HD 37357, etc.) - - the likely evo- 
lutionary precursors to older A stars with debris disks 
like EF Cha - have high signal-to-noise Spitzer spectra 
and numerous solid state features in the 10-25 /im range. 

^ Some of these studies identify solid state features in many 
protoplanetary disks at wavelengths qualitatively similar to phyl- 
losilicate features. However, these authors argue that the f eatures 
are d ue to allotropes of silica (e.g. Sargent et al. 2009; Juhas z~et al.l 
l20lOf) , which we have been able to rule out for EF Cha. 

Their modeling formalism differs from ours in that they adopt 
a "distribution of hollow spheres" to compute dust opacities. We 
use laboratory emission spectra from fine powders. However, their 
computed emissivities quantitatively capture basic features of those 
for powders, specifically a strong, narrow peak at 10 /im. 
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While many (most?) of their disks can be fit by typi- 
cal olivine+pyroxene+silica porous Mie sp ectral models 
and P AH empirical residual templates (e.g. lJuhasz et al.l 
l2010f ). they also provide an interesting future sample to 
apply our spectral modeling to search for phyllosilicate- 
rich protoplanetary disks. 

5. CONCLUSION 

This paper presents analysis of Spitzer IRS data for the 
debris disk surrounding the 10 Myr-old star EF Cha. To 
investigate the composition of the EF Cha debris disk, 
location of the debris dust with respect to the star, and 
size of the emitting grains, we have compared our disk 
spectra to that for disks with a wide range of evolutionary 
stages and analyze it using a sophisticated mineralogical 
model. Our study yields the following main results: 

• Using simple energy balance considerations, the 
EF Cha disk must be a gas-poor debris disk with 
LTE dust temperatures between 225 K and 430 K 
and dust physically located at ~ 1-7.5 AU. Be- 
cause these temperatures are characteristic of equi- 
librium dust temperatures between 0.45 and 1.55 
AU of the Sun in our solar system, we expect that 
the EF Cha debris disk contains terrestrial zone 
material. 

• Unlike most debris disks, the EF Cha disk shows 
evidence for numerous solid state features diag- 
nostic of the composition of parent bodies pro- 
ducing its debris, implying consistent presence of 
small, /iTO-sized grains around the star. Compar- 
ing the EF Cha spectrum with that for other disks 
exhibiting solid state features (mostly protoplane- 
tary disks), shows that the EF Cha disk lacks ev- 
idence for PAH emission and has relatively little 
crystalline olivine and pyroxene. While the general 
shape of the EF Cha spectrum is consistent with a 
high abundance of amorphous silicates, numerous 
maxima secondary features show that additional 
species must be abundantly present. 
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• Minerological modeling shows that the EF Cha de- 
bris disk includes an exceptionally high mass frac- 
tion of phyllosilicates, produced by aqueous alter- 
ation of anhydrous rock, making it unique amongst 
all debris disks studied thus far. Its elevated olivine 
to pyroxene ratio is further evidence of large-scale 
aqueous alteration. The dust may also contains 
water ice and metal sulfides. 



• The amount of dust present is the equivalent of a 
large, ~ 110 km radius asteroid, implying the par- 
ticipation of a large parent body or bodies sourcing 
the observed dust; this sourcing could have been 
caused by collisions in an asteroid belt. 

• Because only a narrow range of coUisional en- 
ergieswill yield detectable debris but not convert 
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tion, EF Cha's debris disk represents a potentially 
rare laboratory for studying phyllosilicate species - 
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Composition of the Best-Fit Model to the Combined Spitzer IRS EF Cha Spectrum 
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TABLE 2 

Composition of the Best-Fit Model to the Spitzer IRS EF Cha Spectrum for Different 

\on POHITIOX^ 
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nod position. With the exception of orthoenstatite, our identification of species from 
not affected by the noise characteristics of the spectra. 


first nod position and then the second 
modeling the combined spectrum are 
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Fig. 1. (Top Panel) Spitzer IRS low resolution (black line) and high resolution (grey line) spectrum of EF Cha. (Bottom panel) 
Signal-to-noise of the low resolution and high-resolution spectra vs. wavelength. The plotted signa-to-noise is smoothed by five pixels in 
wavelength. 
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Fig. 2. — (Left) Extracted low-res spectra for each of tlie two nod positions (black and grey lines). (Right) Histogram plot of the flux of 
the first nod position divided by the flux difference between the nod positions. The distributions are separated by wavelength: 6-15 /im 
(black line), 15-30 /im (grey line), and 30-38 /im (dotted line). 
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Fig. 3. — Comparison between the EF Cha IRS spectral extraction produced from SMART and that produced from the FEPS pipeline 
(top panel). (Bottom panel) The same comparison with the FEPS spectrum rescaled by 0.95, a 5% decrease. 
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Fig. 4. — EF Cha SED including photometric data from the Tycho-II survey and 2MASS (dots) and IRS data (soHd black line). 
The dashed lines correspond to the star + siuglc-tcmpcrature blackbodies; the dotted line corresponds to the combined flux from these 
blackbodies+the star's flux. For the star's flux, we fit the SED to a Kurucz stellar atmosphere model for a Te ~ 7400 K star (grey line). 
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Fig. 5. — (Top) Spectra of a protoplanetary disk (HD 100546), Comet Hale-Bopp, the eje cta from Co met Tcmpcl-1, and several high- 
luminosity debris disks - HD 113766, HD 172555, and HD 69830 compared to that for EF Cha | |Lisse et al . 2008, 2009, 2007b;). Wavelengths 
for PAH emission and other solid state features are identified by vertical dashed lines. (Bottom) Spectra redisplayed in emissivity space to 
better contrast the feature to continuum flux ratios for each object. 
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Spitzer EFCho Disk Spectrol Model 




5 10 15 20 25 30 35 

Wavelength (um) 

Fig. 6. — Spitzer IRS emissivity spectrum of EF Cha with the best-fit spectral decomposition. The central source's photospheric 
contribution has been removed using a Kurucz model with a 7400 K photospheric temperature. Error bars are it 2a. The amplitude of 
each colored curve denotes the relative amount of that species present in the best-fit model (Table 1). For species with no statistically 
detectable emission, the curve is a flat horizontal line. Black points: Spitzer dust excess spectrum, divided by a 600 K blackbody. Orange 
dashed line: best-fit model spectrum. Colored curves: Purples - amorphous silicates of pyroxene or olivine composition. Light blues - 
crystalline pyroxenes: ferrosilite, diopside, and orthoenstatite. Dark blues - crystalline olivine forsterites. Red - amorphous carbon. Deep 
orange - water ice. Light orange - water gas. Yellow lines- phyllosilicates. Olive green - ferromagnesian sulfide (Feo.gMgo.iS). The peak at 
^ 10 fim is caused by a combination of amorphous silicates and phyllosilicates; phyllosilicate species are responsible for structure at ~ 15 
/im, which is more readily apparent in the following flgure. Spectral structure at 16-25 fim is due mostly to forsterite. 
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Spitzer EFCha Disk Spectral Model 
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Fig. 7. — Residual emissivity of the circumstellar dust, after the strong emission due to sihcates has been fit and removed. The remnant is 
dominated by emission from phyllosilicates, amorphous carbon, water ice, and metal sulfides. While there are interesting hints of potential 
water gas emission at 6 fim, and of carbon dioxide gas emission at 15 fim, neither of these features is statistically significant. Emission 
peaks at 9.8 /im and 15 /im due to phyllosilicate species are visible. 
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Spitzer EFCho Disk Spectrol Model 
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Fig. 8. — Same as Figure [6l except with phyllosilicates removed from the model mix of materials and the Spitzer data re-fitted. 
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Fig. 10. — (Top) Olivine vs. pyroxene mass fraction for EF Cha compared to that for other debris disks, protoplanetary disks, comets and 
asteroids. The distribution of ratios define a system evolutionary sequence from unprocessed, primitive rock to highly processed, aqueously 
altered, and pyrolised rock. (Bottom) Atomic abundances relative to solar for EF Cha and other debris disks, protoplanetary disks, and 
comets. 



